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We study the effect of electron-electron interaction on the surface resistivity of three-dimensional 
(3D) topological insulators. In the absence of umklapp scattering, the existence of the Fermi-liquid 
(T^) term in resistivity of a two-dimensional (2D) metal depends on the Fermi surface geometry, 
in particular, on whether it is convex or concave. On doping, the Fermi surface of 2D metallic 
surface states in 3D topological insulators of the Bi2Te3 family changes its shape from convex to 
concave due to hexagonal warping, while still being too small to allow for umklapp scattering. We 
show that the term in the resistivity is present only in the concave regime and demonstrate 
that the resistivity obeys a universal scaling form valid for an arbitrary 2D Fermi surface near a 
convex/concave transition. 

PACS numbers: 72.10.-d,73.20.-r 



Topological insulators (TI) are characterized by a 
gapped bulk spectrum with conducting surface states ex- 
tending across the entire gap. The surface states contain 
an odd number of Dirac cones and are protected against 
any perturbation that preserves time-reversal symmetry 
[T]. A wide variety of interesting physics resulting from 
these surface states is expected to be observed ranging 
from Majorana fermions [2] to magnetic monopoles [Sj. 
Although photoemission and tunneling microscopy [1] 
have convincingly established the presence of such sur- 
face states in these materials, signatures of these states 
in transport measurements are more difficult to observe, 
mainly because of strong conduction in the bulk [5] . With 
recent experimental progress, however, in the ability to 
tune the number of surface charge carriers [5], it is now 
possible to see more clearly evidence of surface trans- 
port. In light of this progress, it is timely to ask what 
is the effect of the electron-electron (e-e) interaction on 
surface transport. Indeed, including e-e interaction is 
crucial for explaining the observed field and tempera- 
ture dependences in quantum magnetotransport [7]. In 
this Letter, we address the manifestation of the e-e in- 
teraction in semiclassical transport within a model of a 
two-dimensional Fermi liquid relevant for surface states 
doped away from the Dirac point. 

An archetypal signature of the Fermi- liquid behavior in 
metals is the dependence of the resistivity (p). With 
an exception of compensated semi- metals [5] , this depen- 
dence in clean conductors arises due to a special type 
of scattering processes- "umklapps" [9l[l0]- in which the 
total momentum of an electron pair is changed by an in- 
teger multiple of the reciprocal lattice vector. Umklapps 
are possible if certain conditions are met, namely, if the 
Fermi surface (FS) is large enough (the band is more 
than quarter full) and if the matrix element of the inter- 
action has sufficient weight at large momentum transfers. 
Otherwise, the umklapp contribution to the resistivity is 



suppressed. In this case, the contribution to p may 
still occur due to the combined effect of the momentum- 
conserving ("normal") interaction among electrons on a 
lattice and electron-impurity (e-i) scattering. Whether 
this really happens, turns out to depend crucially on 
the dimensionality. While the term is allowed for an 
anisotropic FS with a non-parabolic spectrum in three 
dimensions (3D), the conditions in two dimensions (2D) 
are much more stringent In particular, a term 

occurs in 2D only if the FS is either concave or multiply- 
connected; otherwise, the leading e-e contribution scales 
as |T2]. Likewise, the frequency dependence of the 
ac resistivity scales as fl* instead of fl^ [13]. The rea- 
son for such a behavior is that the (17^) term arises 
from electrons confined to the FS contour. For a con- 
vex and singly-connected contour, the momentum and 
energy conservations are similar to the ID case, where 
no relaxation is possible. 

We propose the surface state of a 3D TI as a testing 
ground for the theoretical results outlined above. Pho- 
toemission shows that the surface states of the Se, Te 
based compounds (Bi2Te3, Bi2Se3, and Sb2Te3) have a 
small, singly-connected FS at the center of the BrilUouin 
zone (BZ). Following Fu [T3|, the electronic dispersion in 
these systems can be described by 

e± = ±v/w2fc2 + A2fc6cos2(36l), (1) 

where 6 is the azimuthal angle, v is the Dirac velocity, 
and A is a constant. Corresponding isoenergetic contours 
are presented in Fig. [l] As the Fermi energy increases, 
the FS changes rapidly from a circle to a hexagon and 
then to a hexagram. At some critical value of the Fermi 
energy ep — £c {— 0.16 eV for Bi2Te3, for example [H]), 
the shape changes from convex to concave. Theory [T2 
I13j predicts, therefore, that the e-e contribution to the 
resistivity scales as max{T^, Q*} on the convex side and 
as max{T^, il^} on the concave side. The main result of 
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this Letter is that, near the convex/concave transition, 
the resistivity obeys a universal scahng form 

P = Po+a(-\ e{A)T^ + BT\ (2) 

where pa is the residual resistivity, A = ep — ec, 0{x) 
is the step function, and A, B are material-dependent 
parameters. 

The exponents of 2, 4, and 9/2 in Eq. ([2| are uni- 
versal, i.e., they are the same for an arbitrary 2D Fermi 
surface with a non-quadratic energy spectrum near a con- 
vex/concave transition. We emphasize, however, that 
the surface states of 3D TIs present a unique case of a 
small yet strongly warped 2D FS, where the predicted 
effects can be seen best. The drawback of 3D TIs is 
that they have a large background dielectric constant 
(e 29 — 85 [T5]), and hence, the electron-phonon (e-ph) 
interaction is expected to dominate the T-dependence of 
the resistivity down to very low T[16 . This drawback 
can be circumvented by measuring the frequency depen- 
dence of the optical conductivity at frequencies above 
the Bloch-Gruneisen frequency, where the e-e contribu- 
tion dominates over the e-ph one |17| . 

As in Ref. [12], we adopt an approach based on the 
semiclassical Boltzmann equation (BE) 18 and neglect 
quantum corrections to the conductivity. For simplicity, 
the e-i interaction is accounted for within the l/r ap- 
proximation. First, we consider the dc case (Q = 0). For 
low enough T, when T^e ^ Tei, we solve the BE to lead- 
ing order in the e-e interaction and obtain the correction 
to the residual conductivity as 

^ - ^ / (0 / / / / / t!f V 

X |Mk,p(q)P(Av,)'n(ek)ri(ep)[l - n(ek - u)] (3) 
X [1 - n(ep + w)](5(ek - ek-q - w)5(ep - Cp+q + w). 

Here, t; is the mean free time due to impurity scatter- 
ing, q and lo are the momentum and energy transfers, 
Av = Vk + Vp — Vk-q — Vp_|_q, da\ is the FS element, 
n{£\) — (exp (ei/T) -I- 1)^^, and Mk,p(q) is the matrix 
element of the e-e interaction. To obtain the lowest in 
T term in 6a jj, we project electrons onto the FS, which 
amounts to neglecting w in the arguments of the S func- 
tions. Since the typical values of ek,p and oj are of order 
T, while the typical values of other variables are T inde- 
pendent, Sajj scales as T^, which is the expected Fermi- 
liquid behavior. However, whether the prefactor of the 
term is non-zero depends on whether Av is non-zero 
for all k, p, and q satisfying energy conservation. On 
relabeling p to — p and invoking time-reversal symme- 
try, the arguments of both 6 functions become the same. 
The problem of finding the allowed initial states k and 
p at fixed q now reduces to finding the solutions of the 
equation ek = Ck-q (and the same for k — > p). Geomet- 
rically, this is equivalent to shifting the FS by a vector q 




FIG. 1. (color online). Isoenergetic contours for the spectrum 
in Eq. ([T]). The dashed line corresponds to the critical energy 
for the convex/concave transition. 



and finding k and p as the points where the original and 
shifted FSs intersect. 

Consider first the case ep < Cc, when the FS is con- 
vex. As shown in Fig. [2][a), there are only two points of 
intersection. If k is one of these intersection points then, 
by symmetry, the other point is — k 4- q. Since solu- 
tions for p must belong to the same set {k, — k + q}, the 
scattering process either occurs in the Cooper channel 
({k, — k} — > {k — q, — k + q}) or corresponds to swapping 
of initial momenta ({k, k — q} -> {k — q, k}). In both 
cases, Av — 0, and thus the correction to the con- 
ductivity vanishes. The first nonvanishing term in this 
case is T'*, which can be obtained from Eq. (|3| by ex- 
panding the product of the 6 functions to second order 
in uj. On the other hand, if the FS is concave {ep > Cc), 
there are six possible points of intersection yielding six 
solutions for each k and p [cf. Fig.[2|^b)]. The total set of 
thirty six pairs for {k, p} contains processes other than 
Cooper channel and swapping, and Av is non-zero for 
these processes. 

The analysis presented above is valid either well below 
or well above the convex/concave transition, i.e., when 
|A| ~ ep. We now turn to the vicinity of the transition, 
when |A| ^ ep. First, we focus on the most interesting 
case of A ^ T, when the isoenergetic contours near the 
Fermi energy are concave, and then discuss the case of 
1^1 ^ when both convex and concave contours near 
the FS are thermally populated. Near the transition, sev- 
eral quantities in Eq. ^ exhibit a critical dependence 
on A. First, it is Av which is zero on the convex side 
and non-zero on the concave side. Additionally, there 
are two other quantities which also show a critical be- 
havior. As Figs. [2jc) and [2jd) illustrate, even if the FS 
is concave, it has more than two self-intersection points 
only if it is shifted along one of the special directions and 
the magnitude of the shift is sufficiently small. [These 
special directions are high-symmetry axes that intersect 
the FS at points with positive curvature, as in Fig. 2(b).] 
Therefore, the width of the angular interval near a special 
direction (A0q) and the maximum value of q ((/max also 
depend on A in a critical manner. Approximating J d^q 
by AflqQ^g^^, resolving the S functions, and integrating 
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FIG. 2. (color online), (a) A convex contour lias no more 
than two self-intersection points, (b) For q along a special 
direction, a concave contour has the maximum number (6 for 
a sixfold-symmetric FS) of self-intersection points, (c) Even 
in the concave case, the number of self-intersection points is 
less than the maximal number allowed by symmetry, if q is 
not along a special direction, (d) For q larger than a critical 
value, the number of the self-intersection points is less than 
the maximum number allowed by symmetry. 



over all energies, we obtain 
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X [Av,]L 



^A0q|A/k,,p„( 

Mmaxy I 

kl Pm 1 1 

Vk, • Vp„ • p™ |vk, • q| |Vp_^^ • q| 



,(4) 



where the sum runs over all intersection points, the prime 
denotes a derivative with respect to the azimuthal angle, 
and 1 = 1/|1|. Notice that although a factor of g^ax from 
the phase space of integration cancels with the same fac- 
tor from the 5 functions, it will reappear in the calcula- 
tion of Avj. We are now going to show that 



three roots; the term is non-zero only in this case. 
Since there is a one-to-one correspondence between the 
angles 9* and 6'q, we can find the corresponding interval 
/S.9* instead of A0q. Clearly, the regions on the curve 
where it is non-monotonic are responsible for the multi- 
ple roots [15] • Redefining variables 9 and 9* as measured 
from the invariant points, the non-monotonic part of the 
curve can be conjectured to obey a cubic equation [cf. 
Fig. ic)]: 



a(A)0, 



(6) 



where a(A) cx A and & > is a constant. Indeed, 
we need at least a cubic equation to provide for three 
real roots; whether there is one or three roots depends 
on the sign of a(A) which must be negative/positive 
in the convex/concave regimes, correspondingly. For 
the model spectrum of Eq. we find 6 = 2 and 

«(A) " 



6^(7^)^- The quantity A9* is the ver- 
tical distance between the maximum and minimum of 
this curve which, according to Eq. ^ scales as A^/^. 
Next in line is Avj, which we expand in small q as 



where k2 and ki are any two solutions of the equation 
q • Vk = 0. Referring to the geometry of Fi^. [sl 
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where use of Eq. yields 




we find Av 
Avj cx qA 

Finally, to find gmax, we relax the assumption of small 
q and solve the equation ek-q = eq for arbitrary q. It 
is easier to do this by casting Eq. ^ into an equation 
for the contour in terms of local cartesian coordinates. 
To this effect, we approximate 9* « tan0* ^ dky/dk^ 
and 9 « ta.n9 « —kr^/kp, with kp being the Fermi mo- 
mentum at the invariant point [cf. Fig. |3jd)] and sub- 
stitute into Eq. ^ to get the following contour equa- 
-6fc4/4 + a(A)fc2/2, where kx,y are the mo- 



tion: 



ky — 



menta measured from the invariant points and normal- 
ized by kp . Using this expression to solve for the roots of 



A9„ cx A3/2,, 



cx 



A^/^ and Av^ cx A'^^^ (5) ^k-q = fq one arrives at a cubic equation in k^, which has 



We begin with A9q. Under an assumption (to be justi- 
fied later) of small q, the equation ek — Ck-q = reduces 
to Vk • q = 0, which implies that q is a tangent to the 
FS at the intersection points [cf. Fig. |3]ja)]. Defining 9* 
as an angle between the normal to the FS at any given 
point and q, we plot 9* as a function of the azimuthal an- 
gle 9. Figure |3jb) clearly demonstrates a distinguishing 
feature between the convex and concave contours: 9* {9) 
is monotonic for the former and non-monotonic for the 
latter. The non-monotonic part is centered around cer- 
tain invariant points, i.e., common points for all contours. 
The oscillations reflect the rotational symmetry-sixfold 
in our case-of the FS. From symmetry, if is a solution, 
so is + tt; we thus consider only the domain 9 e [0, tt]. 
We now need to find the angular interval of q about a 
special direction in which the equation 9* = 0q + 7r/2 has 



three distinct real roots if q < 2y/ a{A)/b. This means 
that gmax oc A^/^ (and thus Avj cx A^/^). This, in hind- 
sight, validates the assumption of small q. 

Collecting all the terms together, the resultant energy 
dependence of the critical terms is A0q[Avj]^ cx A^/^. 
which leads to the scaling form of the prefactor of the 
term in Eq. ([2]). The Cooper and swapping channels of 
scattering always contribute a term to the resistivity, 
and the second term in Eq. ([2| accounts for this contri- 
bution. A crossover between the T* and behaviors 
occurs at T ^ epiA/ep)^^'^ ^ cf- A large value of the 
exponent (9/2) indicates that one needs to go sufficiently 
high above the convex/concave transition in order to see 
the term. 

Now we return to the range of energies A < T, when 
both convex and concave contours are populated. Instead 
of a A^/^ factor, we get a factor of T^/^, leading to a 
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FIG. 3. (color online), (a) For small q, points where the 
normal to the FS is perpendicular to q are the points of self- 
intersection (black dots), (b) 6* vs 6 [as defined in panel (a)]. 
Dotted: t < tc\ dashed: e = e^, solid: e > e^. (c) A zoom of 
the non-monotonic part of the graph in panel (b). (d) A por- 
tion of the FS contour showing the geometrical construction 
for the derivation of the equation for the contour. 



r"/2 term in p. This term, however, is subleading to 
the one. Therefore, we conclude that Eq. ^ describes 
the leading T-dependence of the resistivity in all possible 
situations near the transition. 

Next, we discuss the feasibility of observing these pre- 
dictions in an experiment. First of all, one needs to ask 
if the e-ph contribution to the resistivity masks the e-e 
one. In general, the e-ph contribution, which scales as 
at T < Tbc, where Tbg is Bloch-Gruneisen temperature 
(ft! IQK for Bi2Te3 [H]), is expected to be outweighed by 
the (or even T*) one from the e-e interaction. How- 
ever, the e-e coupling may be substantially reduced due 
to high background polarizability of TI materials. In- 
deed, comparing the scattering time of e-ph interaction, 
calculated in Ref. [TB], with that of the e-e interaction, 
we find that the term dominates over the one down 
to a few mK. Even with some uncertainty in the estimate 
of the e-ph time related to screening of this interaction 
by free electrons, the detection of the term in a dc 
measurement seems to be difficult. Instead, as an alter- 
nate route to test our predictions, we suggest measuring 
the optical conductivity as a function of the frequency 
Vt. Indeed, the crossover between the T'^ and forms of 
the dc resistivity is completely analogous to the that be- 
tween the n'' and forms of the optical scattering rate 
T{n) = (a;*)2Re/r?(0) [T3]. It can be readily shown that 
the effective plasma frequency contains the same in- 
tegrals as in Eq. ([s]) [131 |2Q]. Therefore, all the foregoing 
conclusions on the temperature dependence at different 
values oi tp carry over to the frequency dependence. The 
advantage of an optical measurement is that the e-ph part 
of r(ri) saturates [17 for Tbg ^ <C ef, while the e-e 



part continues to grow either as or fi^, depending on 
the sign of A. This advantage was used in the past to 
detect the e-e contribution to T{yi) in noble metals jl7j . 
and we propose to apply the same technique to TIs. 
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